####################################################################################
###
###   Immigration, Public Housing and Support for the French National Front
###   Gloria Gennaro
###
###   Paper Figure A20
###
####################################################################################


library(ggplot2)
library(rdrobust)

################################################################################
# Set Up
################################################################################

# Load Working Directories
source(paste0(wd_main, '/2_code/00_working_directories.R'))

# Load data and clean
source(paste0(wd_code, '/02_data_load_clean_rdd.R'))


data = df_did

################################################################################
# Bdw select
################################################################################

rdbwselect(data$hlm_share_log, data$running_pre, fuzzy=data$sruT, c=0) %>% summary()
df = data[abs(data$running_pre)<898,]

controls = "hlm_share_99 + res_owned_share_99 + imm_share_99 + fn_1995"


################################################################################
# First stage (Figure A20a)
################################################################################

df$res = lm(as.formula(paste0("sruT ~ ", controls)), data=df)$residuals
rdplot(df$res, df$running, p=2, y.label='Policy Application', x.label='Population around the cut-off')
ggsave(paste0(wd_res, '/figures/FigA20a.pdf'), height = 7, width = 5)


################################################################################
# Effect of the policy on housing  (Figure A20b)
################################################################################

df['y'] = lm(as.formula(paste0('hlm_share_log', " ~ ", controls)), data=df)$residuals
rdplot(df$y, df$running_pre, p=2, x.label='Population around the Cut-off', y.label='', title='')
ggsave(paste0(wd_res, '/figures/FigA20b.pdf'), height = 7, width = 5)
